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Using random matrices, we study the reduced dynamics of a two-level system interacting with a 
generic environment. In the weak-coupling limit, the result can be obtained directly from known 
results for purity decay, and result in Markovian dynamics. We then focus on the case of strong 
coupling, when the dynamics is known to be non-Markovian. In this regime, the coupling dominates 
over the local parts of the Hamiltonian, and thus we treat the latter as a perturbation of the former. 

With the help of a linear response approximation, this allows us to obtain an analytical description 
of the reduced dynamics. Finally, we find a transition from non-Markovian to Markovian dynamics 
at a point where the coupling and the local Hamiltonian are comparable in size. 
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I. INTRODUCTION 

In Ref. [1] it was shown that one should expect non- 
Markovian behavior when a central system is coupled 
strongly to a generic environment. In that work every¬ 
thing else but the coupling operator was neglected mm- 
In the present paper, we will study the fate of non- 
Markovianity, when the coupling to the environment is 
still strong, but a local part is also present. The main 
mathematical tool to address these questions is random 
matrix theory (RMT). This theory has found a wide vari¬ 
ety of applications in several fields |1] , including quantum 
chaos, where a direct link between the ensembles stud¬ 
ied in RMT and classically chaotic systems has been well 
established [MS]. Moreover, the idea of complicate inter¬ 
actions, has been extrapolated to encompass interactions 
between two systems, an idea which was formalized, un¬ 
der certain conditions, by Lutz and Weidenmiiller [9]. 
This can be exploited to, say, develop a theory of deco¬ 
herence with RMT; see mmsi. Gonsidering the coupling 
term as the unperturbed system, and the local (free) 
Hamiltonian as the perturbation, we find a critical per¬ 
turbation strength, beyond which the system becomes 
Markovian. At this point the free part is equally impor¬ 
tant as the coupling part. 

While in the infinitely strong coupling case (i.e. with¬ 
out local terms) [T] it was possible to obtain an exact 
analytical solution, here we have to resort to a linear 
response approximation |12j . Even then, the analytical 
solution is quite involved, as it requires the calculation 
of a large number of monomial integrals over the unitary 
group (for simplicity, we will assume the absence of any 
symmetries, including anti-unitary ones ) nanii. 

The paper is organized as follows: In the following sec¬ 
tion, we will describe the system and environment, and 
show that the dynamics of the central two-level system is 
completely determined by a single real function a{t). We 
describe the measure of non-Markovianity which we are 
using, and review known results of the system in the limit 


III 


of strong [T] and weak coupling [IMle]. Next, in Sec. 
we use the results for the evolution of purity to calculate 
the channel for weak coupling. In Sec. EYl we calculate 
the linear response approximation for a{t), when both 
the free part and the coupling term are present. We ob¬ 
tain an explicit expression when the dimension of the 
environment is finite, and a much simpler one in the in¬ 
finite case. We then compare our results to numerical 
simulations. In Sec.jV] we discuss the sharp transition be¬ 
tween non-Markovian and Markovian dynamics halfway 
between strong and weak coupling, in the limit of infi¬ 
nite dimension, where the dimension of the environment 
and the corresponding Heisenberg time are both going to 
infinity. We finish the paper with Sec. lYll in which the 
conclusions are given. 


II. THE SYSTEM 

Gonsider the usual system-environment setting, with 
the Hilbert space being factored in 

H = ( 1 ) 

where TLs corresponds to the system and 'Hf, to the en¬ 
vironment. Moreover, let us choose a single two-level 
system (qubit) as central system, such that dirnHs = 2, 
and a finite dimensional environment with dim He = N. 
The Hamiltonian governing the system is set to be 

H = sl2 G He + V. (2) 

This represents the simplest nontrivial choice for the local 
part of the Hamiltonian, where any dynamics in the qubit 
is neglected. We shall distinguish three regimes: the fully 
coupled system, when s = 0; a strongly coupled regime 
when the norm of the coupling V is comparable to the 
norm of the free Hamiltonian sffej and a weak coupling 
regime when the norm of the coupling is much smaller 
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than that of the free Hamiltonian. The evolution of the 
qubit is given by 


Ps^ = tre 


0 PeU- 


( 3 ) 


where the evolution operator is U* = exp{—iHt). We 
use the Pauli basis, to represent the quantum channel 
induced by Eq. ([^. The corresponding matrix elements 
are given by 

^ ® 0 PeU~*] , (4) 


where (t° = 1 and = ^x,y,z- Notice that choosing 

He and V in Eq. ([^ from unitarily invariant ensembles, 
results in an ensemble of Hamiltonians H that is invari¬ 
ant under local unitary transformations. In the case of 
the central system, this implies that after averaging, the 
channel must be isotropic, so its structure is 


A(‘) = (a(‘)) 


/I 0 0 0 \ 

I 0 a(t) 0 0 I 

0 0 a{t) 0 

Vo 0 0 a{t)J 


( 5 ) 


Here, we introduced the notation (•) for averages over 
the ensemble of random matrices. In the case of the en¬ 
vironment, the above invariance property implies that 
does not depend on the initial state Pe of the envi¬ 
ronment. This allows us to replace pe with the maximally 
mixed state and write 

Q;(t) = ^(tr[cr^ 0 ]leC/“V^ 0 (6) 


A. Full coupling 



FIG. 1. ao{t) as a function of time for different dimensions 
N. Nonmonotonic behavior, causing non-Markovianity, is ob¬ 
served even in the limit of large dimensions, A —>■ oo. Here, as 
well as in all subsequent figures, time is measured in dimen¬ 
sionless units, as explained in the main text below Eq. 


the Haar measure, as explained in min], and obtains 
the general expression 


ao{t) 


4A^|/(f)p-l 
4A2 - 1 


( 9 ) 


where f{t) = ^ exp{—iVjt) is the Fourier transform 
of the spectral density of the Hamiltonian (remember 
that, for s = 0, E is the Hamiltonian of the system). 
Notice that this expression is valid for any unitarily in¬ 
variant ensemble, not just the GUE. One can rewrite the 
above expression as 


4N^l{t) + 2N[l-b2m-^ 

= -7A75-^- 


( 10 ) 


The solution to the fully coupled case, corresponding 
to s = 0 in Eq. ([^, has been worked out in detail in 
Ref. [T]. Here, we only recall the most important results 
as they are to be generalized in the present work. This 
allows us to introduce some notations. For s = 0, the 
quantity to be calculated is 


ao{t) 



0-2 0 1 e 


-iVt 


cr^ 0 1 e 



( 7 ) 


Recall that V is just the coupling term, to be chosen from 
the Gaussian unitary ensemble (GUE) of dimension 2N. 
We shall diagonalize V (and thus the evolution operator) 
with the unitary matrix O. We thus have 

U*=e-*'^‘ = Odiag(e-“^*)OV (8) 


where the {vj}j are the eigenvalues of V. We use units 
for time and energy such that h is eliminated and the 
spectral range of V is equal to (unless stated otherwise, 
the level density for V obeys a semicircle law). As a 
result, energies and times are denoted by dimensionless 
quantities. One then averages with respect to O, with 


where bi is the Fourier transform of the level density of 
V, and 62 is the two-point form factor without unfolding; 
cf. Ref. [IH]. For the GUE, both functions are known 
analytically and are given in Appendix [A] in Eqs. (A3) 
and (A4) (together with further details). 

Spectral correlations are expected to be limited to an 
energy scale of the order of the mean level spacing, which 
is N times smaller than the energy scale of the level den¬ 
sity. As a consequence, the relevant time scales for 61 and 
62 become very different for large dimensions. We chose 
matrices V from the GUE such that (VijVki) = SjkSu/N. 
In this way, in the limit A —> 00 , the level density tends 
to a semi-circle on the interval (—1,1). As we have set 
h = 1, the relevant timescale for 61 is therefore of order 
1 (we call this timescale “macroscopic”), while for 62 the 
relevant timescale is the Heisenberg time which is of or¬ 
der N. In the limit N ^ 00 , we get for the GUE an 
oscillating function in time: 


lim anU) 

AT-S-OO 


t 


( 11 ) 
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B. Non-Markovianity in the fully coupled system 

Quantum non-Markovianity does not have a unique 
definition. Definitions include considering any deviation 
from the Lindblad equation as non-Markovian behav¬ 
ior m, the backflow of information from the environ¬ 
ment into the system |20j . and also the impossibility of 
defining an instantaneous quantum map for intermediate 
times m- Accordingly, several measures have been pro¬ 
posed to quantify the degree of non-Markovianity, each 
with different properties and problems [22l [23]. How¬ 
ever, for simple channels, like a depolarizing channel, as 
in our case, most definitions coincide as far as the dis¬ 
tinction between Markovianity and non-Markovianity is 
concerned [T| , even though the measures of the degree of 
non-Markovianity are usually not comparable. For the 
sake of definitiveness we shall use the measure proposed 
in (20), although other measures could be easily incor¬ 
porated in this framework. The measure is defined as 
the maximum of the integrated backflow of information 
measured in terms of increasing distinguishability, where 
the maximum is taken over all possible pairs of initial 
states. In the present case, where the quantum process 
is determined in terms of the function a{t), one gets [T| 

poo 

M = 2 dt d{t) . (12) 

J d>0 

The measure will be greater than zero if and only if 
d{t) > 0 for some time, i.e. if the Bloch sphere ex¬ 
pands during a time interval. One of the results of 
Ref. [T] says that the system will generically display non- 
Markovianity, even in the limit of an infinite dimensional 
environment [N —^ oo); see Fig. One may compare 
the present model to the case of an environment mod¬ 
eled by a collection of harmonic oscillators, characterized 
by a spectral density J{w)-, see [21|, chap. 10. In those 
models this spectral density has a role similar to the level 
density in ours (however, see jlSlIini)) as it is the forms 
of those functions which determine the reduced dynam¬ 
ics and thereby the (non-)Markovianity. This similarity 
is surprising, since we are dealing with a very strong cou¬ 
pling limit, whereas the description based on the spectral 
density relies on a weak coupling approximation. In this 
respect, we also find it surprising that at strong but finite 
coupling, our model shows a transition to Markovian dy¬ 
namics, independent of the level density. That case will 
be discussed in Sec. IIYI 

At finite N, the non-Markovianity has two contribu¬ 
tions acting at different time scales. The first comes from 
the oscillations in the one-point function 6i(t), which ap¬ 
pear on a timescale independent of the dimension N of 
the environment (the macroscopic timescale). The sec¬ 
ond contribution comes from a recovery of a{t) between 
the first zero of Ji{2t) and the long-time limit 

That occurs on the timescale of the Heisenberg time th 


of the environment, which is proportional to N. In the 
semiclassical limit, N —>■ oo, the Heisenberg time goes to 
infinity and the recovery goes to zero. As a consequence, 
limjv_>.oo M =0, also. 


III. THE WEAK-COUPLING LIMIT 


The behavior of a(t) in the weak-coupling limit can be 
deduced from previous results [HI [16] , where the purity 
for a model equivalent to Eq. © was studied. In that 
limit, the relevant time scale is the Heisenberg time th of 
the Hamiltonian He of the environment. Since the focus 
was then put on the evolution of purity, P = tr , we 
use the fact that purity can be expressed in terms of a(t) 
from Eq. ([^ as follows: 

P(t)=tr(A(*)[p])' = i±|i^. (14) 

Switching from the parameter s, which scales the factor¬ 
ized term, to A scaling the coupling, we can write 

H = He®l + XV. (15) 

V is chosen from a GUE, but now with an iV-independent 
scaling (VijVki) = SjkSu. In order to map this Hamilto¬ 
nian on Eq. ([^, we would have to set A = l/(sA^). In 
the case that He and V are both members of a GUE, 
it was found that, in the linear response approximation, 
the average purity is given by 

-Plr(Q = 1 - X^git) ( 16 ) 


with 

g{t) = 2tmax{t,TH} -b ^(min{t, th})^, (17) 

3rH 

with Th being the Heisenberg time of He- 

To go beyond the reach of linear response theory, we 
exponentiate the result, such that (i) the first two terms 
in a Taylor series coincide with the linear response result 
and (ii) the asymptotic value coincides with a theoretical 
expectation. Such a heuristic procedure, known as expo¬ 
nentiation, has lead to excellent results [H]. In our case, 
the procedure leads to 


^elr(4 = 


1 

2 


2 exp 


liPLKit)-!) 


(18) 


Relying on self-averaging, which is often the case in these 
kind of systems [16], one can reconstruct a{t) for moder¬ 
ate values of the perturbation. Thereby, we obtain 


A 2 

a{t) = exp ( — -g{t) 


(19) 


Notice that when the Heisenberg time becomes infinite, 
we obtain an exponential decay for a(t), a result also 
known as the Fermi golden rule. Notice also that g{t) as 
defined in Eq. (17) is a monotonically increasing function. 


which implies that a(t) is monotonically decreasing. This 
means that the corresponding dynamics is Markovian, 
independent of the shape of the level density. 
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IV. THE STRONGLY COUPLED SYSTEM 


equipped with the Haar measure. Eq. (24) may be rewrit¬ 
ten as 


So far we found that, for iV —>■ oo, the fully coupled 
system (s = 0) shows non-Markovian dynamics, while at 
weak coupling, the system becomes Markovian. In this 
section, we consider the crossover region, when s is small 
but finite. The linear response theory developed below is 
applicable as long as s 1 when the free evolution term 
and the coupling in the Hamiltonian in Eq. are of the 
same size. From a technical point of view, the calculation 
is much more demanding than usual, because the linear 
response expansion is around the fully coupled case. 


A. Linear response theory 
We will calculate 

a{t) = ^ (tr [a^ (g) ]lee”*-^V^ (g) ) (20) 

with the ensemble defined in Eq. (|^. To apply lin¬ 
ear response theory for small s, we consider the unper¬ 
turbed propagator to be and the perturbation siJg. 

Hence, we have for the echo operator: 

e^vt e-Ht « 1 _ js f dt'Heit') 

Jo 

-s^ f dt' [ dt"H,{t')Hein, (21) 
Jo Jo 

where X{t) = e*^*Ae“*^* denotes the interaction pic¬ 
ture of operator X. After some calculations, detailed in 
Appendix [B| we find that 

a{t) ao{t) - s'^a 2 {t) (22) 

where 

a2(i) = J dt" (23) 

and 

(24) 

Ad) = (tr[e*^d-‘")a,e-*^d-Oiy^e-*^‘'a,e*^‘"iIe]). 

(25) 


B. Averaging over the unitary group 

We shall work in the eigenbasis of the environmental 
Hamiltonian, so that He = diagefc. Let us call O the 
matrix of eigenvectors of V so that e*^* = with 

V being the eigenvalues of V. Since V is taken from a 
GUE, O must be an element of the unitary group U (2 V) 




In this equation the Einstein summation convention is 
used. The indices a, b, c, and d run through the basis 
states of the qubit; the indices i, j, k, and I through 
those of the environment; and the greek indices through 
the 2N eigenstates of the coupling operator V. Equa¬ 
tion (26) is composed of three independent parts: The 


first part contains time and the eigenvalues of the cou¬ 
pling. The second one, contains the eigenvalues of the 
environment Hamiltonian, and the third part contains 
the term (—)“+‘^ and the eigenvectors of V. The term 
Ad) can be decomposed similarly. Notice that one can 
go from Eq. (24) to Eq. (25) performing the following 
substitutions: 


Ad) 

t 

t' 

t” 


A(2) 

t" 

t 

t' 


Ad) 

a 

/3 

7 

5 


Ad) 

5 

a 

P 

7 


(27) 


Using these rules, one can w rite the analogous expression 
for Ad)j starting from Eq. (26). As is well known m, 
averages over the unitary matrices with respect to the 
Haar measure are invariant under arbitrary permutations 
of columns and/or rows. Hence, the result of those av¬ 
erages only depends on the question of whether these 
indices coincide among each other or not. We may use 
this invariance property to get rid of the factor (—)“+‘^ as 
follows: Assume i 1; then the row ai is always different 
from dl and the group average does not depend on a and 
d, so the summation over a and d can be factored and 
yields J2a — 0- Therefore, we may restrict the 

summation to the case i = 1. 


The different terms in the summation in Eq. (26) can 


be grouped according to the degeneracy of the indices; 
the particular value of each index is unimportant. One 
can therefore divide the set of values for the four Greek 
indices into 15 different partitions, which will be enumer¬ 
ated as follows: 


1 

a = /3 = 7 = (5 

2 

a = /3 = 7 ^ (5 

3 

a = /I = (5 7^ 7 

4 

a = 7 = (5 7^ /3 

5 

a ^ P = J = 5 

6 

a = /3 7^ 7 = (5 

7 

a = 7 7^ /3 = (5 

8 

a = (5 7^ /3 = 7 


9 

a 

= 

/3 

7^ 

7 


6 

10 

a 

= 

7 

7^ 

P 


6 

11 

a 

= 



/3 


7 

12 

a 

7^ 

6 


/3 

= 

7 

13 

a 

7^ 

7 


/3 

= 

5 

14 

a 

7^ 

/3 


7 

= 

6 

15 

a 

7^ 

/3 


7 


6 


(28) 


For the latin index pairs we can proceed likewise. Due 
to the invariance properties of the averages of the mono¬ 
mials, based on the above labeling of the partitions, we 
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can write 


7^1 J^l 

^( 2 ) = C'^ F^‘^K 


(29) 


Notice that we are using capital latin letters as indices 
for the different partitions. In this equation, C is a vec¬ 
tor containing all Cj cases, in which the terms {sjEk) and 
are taken into account; in the matrices 
the group averages over the monomials of matrix ele¬ 
ments of O are arranged, and the time-dependent phases 
containing the eigenvalues of the c oup ling are are in¬ 
cluded in The partitions, Eq. (28 1 , with respect to 


row indices (latin index pairs) and column indices (greek 
indices) have different multiplicities, which are included 
in the vectors C and respectively. The factors C/ 

are the same for A^^'> and . We find that 

Cl = -C 2 = -Cs = -Ce = -Cr = 27V, 

C 3 = C 4 = -2{N - 2), Cg = 2N{2N - 1) 

Cg = 4(iV - 1), Cio = -C 12 = Ci 3 = Ci 4 = 4(iV - 1) 
Cii = -4(7V - 1){N - 2), Ci 5 = - 1)(7V - 4). 

The group averages appearing in the matrices are 

calculated exactly for arbitrary TV, based on recursion 
formulas developed in [ 1 ^, available as computer code 
in m- We report the results of the vectors 




1 ("v 1 2(TV-l)(TV-t2) TV-1 2(TV2-h 3TV-b 1) TV - 1 

TV(2TV-bl)(2Af-b3) 2TV- l’ 2TV - 1 ’ {2N-1)' (2TV - 1) N{2N-1)' 

{N-1){N+ 2){2N+ 1) TV-1 TV-1 3TV-b2 A^-l 

Af(2TV- 1) ’“TV(2TV- 1)’“2TV(2TV- l)’“2TV(2Af- l)’“2Af(2TV- 1)’ 

TV - 1 4TV3 -b 6TV2 -3TV-2 TV-1 5 \ 

“2TV(2TV- 1)’ 2TV(2TV- 1) ’ “2N(2TV - 1) ’ 2(2Af - 3)(2TV - 1)7 ’ 


and 


1 / ^ TV-1 A^-l TV-1 A^-l A^-l 2(Af-l)(TV-b 1) 

“ N{2N -b 1)(2TV -b 3) V 2N - 1 ’ 2TV - 1 ’ 2TV - 1 ’ 2TV - 1 ’ ~ TV(2TV - 1) ’ Af(2TV - 1) ’ 

(TV-b l)(4N-b 1) 2TV2-b2TV-bl (TV - 1)(TV-b 1) 2TV2-b 2TV-b 1 2TV2-b 2TV-b 1 (TV - 1)(TV-b 1) 
Af(2TV- 1) ’ 2Af(2TV-l) ’ TV(2TV - 1) ’ 2TV(2Af - 1) ’ 2TV(2TV - 1) ’ TV(2TV - 1) ’ 

2TV2-b2TV-bl 2(Af-l)(TV-bl) \ 
2TV(2TV-1) ’ (2TV-3)(2Af-1)7' '' ’ 


C. Average over the eigenvalues of V 


the help of the auxiliary functions 


We now calculate the components of the time- 

dependent factors As we are mainly interested in 

the case of large TV, we shall ignore all spectral correla¬ 
tions, as these could only affect the dynamics of the qubit 
at times proportional to TV, where a{t) already is of order 
of 1/TV. We have seen this explicitly in Sec. IIA where 
we considered the case of full coupling, s = 0. We expect 
that, for the perturbed case with finite s, the situation 
will be similar, and will be justified a posteriori with the 
numerical simulations. In other words, we assume that 
the eigenvalues of the coupling term V have a semicircle 
spectral density, but are otherwise statistically indepen¬ 
dent. Note, however, that in principle, one could take 
into account correlations and describe the behavior up 
to times of the order of the Heisenberg time, if required. 

That said, all components Fj ' ' will depend solely on 
the Fourier transform of the spectral density bi. With 


Fix) 

Qix,y,z) 

nix,y,z) 


i2N)\ 
i2N -2)\ 
i2N)\ 
(2A-3)! 

(2TV)! 

(2A^-4)! 


blix), 


bi{x)bi{y)biiz), 
bi{x)biiy - x)biiz 


(32) 


y)bi{z)- 
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TABLE I. (Color onlin e) Order of magnitude of the 15 differ¬ 
ent terms [see Eq. (28l] contributing to Eq. (22|. The leading 
terms are shaded (in green). Thus, for large A^, it is enough 
to consider the 13th term of and the 15th term of 


we may write 


= 2N, 

1 

II 

= Fin, 

1 

1 

+ 

II 

i-H 0 

(E3W)=E(t'-n, 

(EW)=g(t-t",t'-t,t"-0, 

(EW)=E(t-t'), 

{Fil'>) = git,t" -t,n, 

(E^')) = Fit), 

=g)t,-t + t' -t' 

II 

{Fil^) = git,t' 

{F^^^)=Fit-t'+n, 

{F^l'>} = nit,t',n- 




(33) 


Finally, using the mapping (271, we also obtain the com¬ 
ponents of : 


(i^f)) = 2iV, 

)) = F{t"), 

{Ft^)=Fit"-t), 
{Fi^'>)=F{t"-t'), 

-t + t'), 

) = F{t), 


{F[f)=g{t'\t-t\t), 

{F[f) = g{t" + t' 

{F^f) = g{t" -t), 

{F[f) = n{t",t,t'). 

(34) 


Equations (22), (23), and (29), together with Eqs. (30) 
to (34), provide the final, general result. It is valid, either 
in the absence of spectral correlations in iLe, or for large 
N at times sufficiently small compared to the Heisenberg 
time. In our case, bi{t) is given in Eq. (11), which cor¬ 
responds to a semicircle level density. However, other 
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FIG. 2. (Color online) Exact value of the integral 02 ( 1 ) in 
Eq. (231, with all terms (yellow curves) and only taking into 
account the leading terms in N (blue curves); that is, the 
value of for several values of the dimension N. 


cases with different level densities could be considered, 
also. Our general result still requires the evaluation of 
the double time integral in Eq. (23). Typically, one would 
have to do this evaluation numerically. 


D. The solution for large dimensions and times 


It is possible to simplify the general expressions, dis¬ 
cussed above, considering explicitly the limit of large N. 
Table |l] indicates the leading order in N~ ^ of all relevant 
terms, in Eqs. (30), 111 ), and ( [^ . By proper selec¬ 
tion of the highest order terms, we obtain for a^(t) = 
limjv-foo c( 2 {t) the following 


a“(t) = 2 ^ dt' -t') 

-b,it")b,it-t")bi{t'-t)b,it')]. (35) 


We have tested the reach of this limit in fig. where 
we can see that already for N = 2^° there is almost no 
visible difference, for the times reported, between the full 
expression and the large-environment limit. Although 
this expression means a considerable simplification for 
the 61 from a semicircle level density, we were still not 
able to solve the time integrals in closed form. We found 
only one case where that is possible. This case, where 
the level density is a Gaussian function, is treated in 
Appendix 


E. The solution for large times 


Linear response theory is valid whenever the correc¬ 
tions in the echo operator, Eq. (21), with respect to the 


identity are small; that is, whenever |e*^‘e — 1 | <C 1 
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FIG. 3. (Color online) Comparison of the linear response 
theory, with an exponential tail, and the numerical simula¬ 
tion, with N = 64, for several perturbations. 


(here, | • | denotes the operator norm). This implies that 
the eigenvalues of the echo operator must remain close 
to 1. Departure from that happens for any value of the 
perturbation s, for sufficiently long times. However, the 
smaller the s, the larger the time of the validity of the 
linear response approximation. 

The extension of linear response formulas in this con¬ 
text is often done using exponentiation, as in Sec. m 
However, in the present case such attempts have been 
unsuccessful [5S]. As an alternative, we have opted to 
combine the two linear response theories, namely the ones 
discussed in Sec. IHII and Sec. IIV Al We shall use the lin¬ 
ear response formula Eq. (221 until the time in which the 
largest (in absolute value) eigenphase reaches ±7r. After¬ 
wards, an exponential decay is fitted. 


V. THE TRANSITION FROM 
NON-MARKOVIAN TO MARKOVIAN 
BEHAVIOR 


As the coupling of the system diminishes (that is, when 
s increases), one should fall back to the Markovian case in 
the large-dimension limit [mm]; cf. also Sec. dTO This 
is indeed observed in fig. where the curves for a{t) 
become monotonic as s increases. Thus, we would like 
to know whether there is a particular value for s beyond 
which the dynamics is Markovian. This question is an¬ 


swered in fig. 


^ where the measure for non-Markovianity, 


from Eq. (121 is plotted as a function of the coupling s. 


The points mark the numerical results for A4 where the 
integration in Eq. (12) has been restricted to the interval 
t G [0,10]. While this introduces an error at small di¬ 
mensions, this error vanishes at large N. The solid lines 
mark the same quantity, but calculated from the com¬ 
posite linear response results shown in fig. One can 
observe that there is a seemingly sharp transition in the 
large-dimension limit, which is not observed for smaller 
dimensions due to the two-point correlations that cause 
an increase in the function a, and are not taken into 



FIG. 4. (Color online) Measure of non-Markovianity for a 
random Hamiltonian of the form Eq. §. Here, both matrices 
V and He are taken from the CUE, and the dimension of the 
environment N is varied. Symbols indicate calculations done 
with the linear response theory, extended with an exponential 
decay, while continuous lines are obtained numerically. One 
can see that, at a critical intensity of the coupling s ~ 0.4, 
the system switches from non-Markovian to Markovian. It 
should also be noted that, as the dimension increases the ap¬ 
proximations are more accurate. 


account in the linear response results. 

It is remarkable that a critical value of the coupling 
is needed to go from one regime into the other. It must 
be noticed, however, that in our calculation we are using 
an ensemble of Hamiltonians to describe the environment 
and the coupling to it. In a real experiment this would 
correspond to measurements which require many repeti¬ 
tions of the quantum process, during which the dynamics 
in the environment changes, e.g., due to fluctuating clas¬ 
sical fields. A particular member of the ensemble will 
exhibit oscillations that will result in non-Markovianity. 
However, one should distinguish oscillations due to the 
particularities of the system, from generic oscillations due 
to general features of the whole ensemble. 

For completeness, we have also studied the case in 
which the qubit has an internal Hamiltonian, where 
Eq. (§ is substituted by 

H = ujaz + st2 ® He + V. (36) 

This Hamiltonian is no longer invariant under unitary 
transformations in the central system, and hence Eq. ([^ 
is no longer valid. Instead, the new quantum channel will 
be a combination of a dephasing and a depolarizing chan¬ 
nel. The only energy scale retained in the limit of large 
dimensions is the spectral span of the coupling Hamilto¬ 
nian V [the level density has the shape of a semi-circle 
in the interval (—1,1)]. Therefore, one may expect that 
the effect of the additional term depends on the size of 
w as compared to the spectral span. Hence, for w ^ 1 
the effect should be negligible, we do expect differences 
for w > 1. In Fig. we present simulations for differ¬ 
ent values of w. The figure shows our measure for non- 
Markovianity as a function of s, just as in Fig. 4. We can 
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FIG. 5. The effect of an internal Hamiltonian in the central 
system, as in Eg. |36[ is studied. The level splitting, w, ampli¬ 
fies the non-Markovian effects, but apparently conserves the 
transition from non-Markovian to Markovian behavior, which 
occurs at approximately the same critical value for s. 


observe, that increasing uj leads to larger values for the 
measure, but that the transition from non-Markovian to 
Markovian behavior is essentially unchanged. 


VI. CONCLUSIONS 


We considered a quantum two-level system coupled 
to a generic environment modeled by random matrix 
theory. We obtained analytical expressions for the re¬ 
duced dynamics using linear response approximations, 
both around the weak and strong coupling limits. For 
the weak-coupling limit, an explicit expression is ob¬ 
tained. The corresponding expression involves a double 
time-integral of one- and two-point functions of the cou¬ 
pling in the general case, a of one-point functions in the 
large dimension limit. In the limit iV —>■ oo, the result 
becomes much simpler [see Eq. (35)]: only two terms 
survive, which contain one-point functions. Neverthe¬ 
less, the double time-integral still has to be evaluated 
numerically. 


We then studied the degree of non-Markovianity in 
the system, using the measure proposed in [5D], based 
on distinguishability. We show that the degree of non- 
Markovianity of the s = 0 case (infinite coupling) con¬ 
sidered in Ref. [1] diminishes as the coupling term be¬ 
comes less important, and that in the large-V limit it 
vanishes at a point where the free and coupling terms of 
the Hamiltonian are of equal size. 


Acknowledgements- Support by projects CONA- 
CyT 153190, CONACyT 129309 and UNAM-PAPIIT 
IN111015 are acknowledged. 


Appendix A: Details for the fully coupled case, in 
the GUE case 

The spectral correlations for the GUE are expressed in 
terms of the functions 


-2NE'^/4: 

= =n,{EVN), 

^21j\^2e/{2N) 

where T-Lj are Hermite polynomials [18) . 

We fix the normalization so the average of the square 
of the diagonal elements in the matrices is 1/N. Then, 
the spectral density for finite dimensions is 

2N-1 

R^{E) = Y, HE? (Al) 

i=o 

and the cluster function, containing the correlations be¬ 
tween levels, is 


'2N-1 

T2{E,,E2) = I 5] UE,^,{E2) 

j=0 


(A2) 


If we define 


blit) = ^ / dEe-^^^R,{E) (A3) 

1 

b2it) = ^J dE,dE2e-d^^-^-^^T2iEi,E2) (A4) 


we have, for this case, 


{aoit)GVE) = , (A5) 


Since 


N^fit)\^)=N+ 

j AEAE2e~"^^^~^^'^\RiiEi)RiiE2) - T^iE^.E^)]. 


In the large dimension limit, we have 

Jim 


blit) = 


t 


(A6) 


(A7) 


and thus 


lim (ao(t))GUE = 

N—^ao 


Ji( 2 t) 

t 


(AS) 
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Appendix B: Details for the linear response theory 


Now, we write a{t) in terms of the echo operator as 
follows: 


^{t) = tr 


(72 (g) 1 e * M CT 2 (g) 1 e*^ * 

Using Born approximation Eq. ( |21[ ), we obtain 
a{t) = ti[az{t) M (Tz (g) 1 ] 

« ao(0 ~ f dt"tTA{t",t',t), 

Jo Jo 


(Bl) 


where 


ao{t) = tr [dzit)az g) 1] (B2) 

represent the exact known solution for s = 0, and 
A{t”,t',t) = dz(t)Heit')He{t”)az - az{t)He{t')azHe{t'') 

+ - dzit)He{t'')azHeit'). 

Due to the fact that the matrices Hg, dz and ctz g) 1 are 
Hermitian, we obtain the useful identities 

tY[dz{t)He{t')He{t")az\* = ir[azHe{t")He{t')dz{t)\, 


a variance equal to their inverse dimension. Let H be 
either V or Hg, so the normalization condition reads 

{H,, Hu) = ^ . (Cl) 

That implies for the average of the trace of and for 
the square of the trace of H 

= ^=N, (C2) 

jk 

[triH)f = J2{H,,Hkk) = = 1- (C3) 

jk j 

In turn, this means that the eigenvalues Sj of He lie es¬ 
sentially in the interval (—2,2), and have a semi-circle 
distribution, for large N. In addition, since 

(ef) = 1, (C4) 

and, since J^jki^j^k) = N + N{N - l)(ej£fc) = 1, 

{£jSk)j^k = (C5) 


Appendix D: The Gaussian PUE 


and 

tT:[dz{t)He{t')azHe{t'')\* = tv[He{t'')azHe{t')az{t)] 

= tY[dz{t)He{t")azHe{t')] ■ 

This implies that the trace of A{t”,t',t) can be written 
as twice the real part of the trace of Ac(t",t',t), where 
the latter quantity only contains the two terms on the 
left-hand side of the above equation. We may thus write 
for a(t) in the linear response approximation: 

a(t) « ao(t) — 2 Re f dt' f dt" tr[ Ac(t", t', t) ] , 
Jo Jo 

where 

Ac{t'',t',t) = az{t)He{t')He{t")az-dz{t)Heit')(TzHe{t"). 
Now we split Ac in its two parts, 

trAc(t",t',t) = - AC), 

insert identity operators e“*'^*e®'^* wherever necessary, 
and use the cyclical property of the trace, to rewrite more 
conveniently the term under the integral. 


Appendix C: Normalization of the ensembles 
considered 


In Sec. IV A we are free to consider any normalization 
condition. 

We conveniently assume that, with respect to an arbi¬ 
trary basis, the matrix elements of both V and He have 


Initial calculations were done in a Poissonian ensem¬ 
ble with Gaussian level density. This could correspond 
to spin models. Even though non-Markovian effects are 
not visible here (the Fourier transform of a Gaussian, 
is another Gaussian), some results are easier to obtain, 
and provide a guide to the more complicated calculations 
in the GUE case. We present some details here, which 
might also provide a guide for the calculations using other 
spectral densities. 

We now assume that the coupling V is taken from the 
GPUE (Gaussian PUE). This means that V is chosen as 

Ugpue = UDU^ , (Dl) 


with U chosen from the unitary group with the Haar 
measure, and D is a diagonal matrix with Gaussian in¬ 
dependent numbers with a — 1. This means that the 
variables in Eq. (26) are, again, independent Gaussian 


variables with mean zero and unit standard deviation. 

The subseq uen t calculation runs i n id entical way as 
shown in Sec. IV except that in Eq. (A71, we have 


bi{t) = e” 


(D2) 


and one has to propagate this difference through Eq. ( [^ . 
In this case, some of the integrals involved in the terms 
can be performed, though not all. 

In this case, the compact expression for the large di¬ 
mension limit. 


aTit) 


V^te Erf ( ^ 


7re 2 Erf 



, (D3) 


is obtained. 
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